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The order and the nature of the finite-temperature phase transition of QCD with two flavors of dynamical 
quarks is investigated. An analysis of the critical exponent of the specific heat is performed through finite-size 
and finite-mass scaling of various susceptibilities. Dual superconductivity of QCD vacuum is investigated using a 
disorder parameter, namely the v.e.v. of a monopole creation operator. Hybrid R simulations were run at lattice 
spatial sizes of 12^, 16^, 20^ and 32"^ and temporal size Nt = 4, with quark masses in the range arriq = 0.3 — 0.01. 



1. Introduction 

In QCD with two dynamical fermions a fi- 
nite temperature phase transition is expected 
at iriq = 0. The basic question of the or- 
der of this phase transition has been exten- 
sively studied in the continuum by an effec- 
tive sigma model analysispP and on the lat- 
tice with Kogut-Susskind pi3l4l5l6| and WilsonU] 
fermions. Nonetheless a conclusive answer has 
not yet been found. 

At zero quark mass in the Pisarski-Wilczek sce- 
nario, if the axial anomaly is effectively restored 
the transition is predicted to be a fluctuation- 
driven first order transtion; in the opposite case 
the transition can be second order in 0(4) uni- 
versality class; mean field critical behavior with 
logarithmic corrections is also possible as a limit- 
ing case between the previous two. Since first or- 
der phase transitions are generally robust against 
small perturbations, we expect the first order be- 
havior also if the effective coupling of the anoma- 
lous current is small enough. K-S fermions break 
explicity the 0(4) simmetry to 0(2) at finite lat- 
tice spacing, so in the case of a second order phase 
transition we expect 0(2) critical behavior until 
the continuum limit is recovered. 

The mass term in the QCD Lagrangian corre- 



sponds to an external magnetic field in the effec- 
tive sigma model, thus we expect the transition to 
remain first order at finite and small quark masses 
if the transition were first order at rrig = 0. In the 
other cases, the mass term completely washes out 
the phase transition to an analytical crossover. 

To determine which of the above possibilities is 
actually realized, the standard test of universality 
is a comparison of the critical exponents. Critical 
exponents can be extracted from the scaling anal- 
ysis of susceptibilities obtained from the second 
derivatives of the partition function: 
d 
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L^Nt. Do is temporal component of 
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Xm = 

Xij — 
Xee = 
Xie = 

with V = ij^i^t, ^0 
the Dirac operator and P^, Pt indicate the aver- 
age spacial and temporal plaquette respectively. 
Cv denotes the specific heat and it can be ex- 
pressed in term of Xij, Xee and Xie', Xm is the 
susceptibility of the chiral condensate and only 
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Figure 1. p as a function of L for Runl. Reweight- 
ing curves are also shown together with MC data. 

its disconnected component is taken into account 
in the present work. The first of these suscep- 
tibihties p is defined in term of the disorder pa- 
rameter introduced in |8I9| . jj, creates a magnetic 
charge defined by some given abehan projection 
and (//) detects magnetic monopoles condensa- 
tion. The parameter has been tested in quenched 
theoryPEI and it is weU defined in full QCDfTU]. 
The independence from the chosen abelian pro- 
jection has also been successfully tested |^ITT] . 

The peaks of these susceptibilities define the 
pseudocritical coupling and their heights obey 
scaling laws. Since we have two parameters, 
namely the bare quark mass am,q and the finite 
spatial size L, in order to obtain scaling relations 
of just one variable, we constrained our parame- 
ters by the relation am^-L^'' = cxjst. In proximity 
of the critical point, the singular part of the free 
energy density fs satisfies 

fsit,mg,L-^)^b-''UbyH,by'^mg,bL-^) (7) 

where t = [T — Tc)/Tc is the reduced temper- 
ature, yt and yh are the thermal and magnetic 

Table 1 

Run parameters for the numerical simulations. 



Ls ■ TO,r varies in the range 


8.9-15.8] 




aniq 


# Traj. 


Ls 


Runl 


Run2 


Runl 


Run2 


12 


0.153518 


0.307036 


22500 


25000 


16 


0.075 


0.15 


38000 


54700 


20 


0.04303 


0.08606 


12800 


12500 


32 


0.01335 


0.0267 


7000 


15100 



critical exponents, and b is an arbitrary scale fac- 
tor. From this and using the above constraint, we 
obtain the scaling laws 

Pc{L-^) = PM + cpLy* (8) 
p'it^L) = Ly'^itLy) (9) 
C{?"^ - co + ciL^^'-'^ (10) 

xUt.L) = L^y^-'^tLy-) . (11) 

2. Numerical Results 

We used the plaquette action for gluons with 
the Kogut-Susskind quark action. Configura- 
tion updating was performed adopting the stan- 
dard Hybrid R algorithm. The lattice tempo- 
ral size was fixed at A''^ = 4, the molecular dy- 
namics step size was Ar = 0.005 and the molec- 
ular dynamics trajectories were generated with 
unit lenght. As discussed in the previous sec- 
tion, we constrained simulation parameters by 
aniq ■ Ly*^ = cost. Expecting 0(4) scaling from 
previous works, we fixed the magnetic critical ex- 
ponent to be yh — 2.49; notice that the 0(2) value 
is also very near to that. We performed two paral- 
lel set of runs, referred to as Runl and Run2, the 
only difference between the two beeing that for 
Runl we set aruq ■ = 74.7 whereas for Run2 
aniq ■ = 149.4. Tabled summarizes the sim- 
ulations parameters for our runs and shows the 
statistics collected from each run. Results from 
Runl suffer from the low statistics collected at 
L = 32 and should be considered as preliminary. 

The behavior of the p susceptibility is shown in 
Fig-E For T < Tc p is size indipendent showing 
that (/i) 7^ 0; for T > Tc p shows a growing linear 
with L showing that (/i) is strictly zero in the 




Figure 2. Pseudocritical coupling vs 1/Ls 
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Table 2 



yt 


(P) 


yt (x<t<t) 


yt (A) 


yh (Xm) 


Runl 


Run2 


Runl 


Run2 


Runl 


Run2 


Runl 


Run2 


3.24 ±0.05 


3.15 ±0.03 


2.13 ±2.48 


2.37 ±0.69 


1.51 ±0.14 


1.35 ±0.1 


2.78 ± 1.27 


3.17 ±0.48 



limit i — > cx). 

A power law fit of the p peak as a function of L 
gives a value for yt (Eq. O . Other estimates for 
this exponent can be obtained from analogous fits 
(Fig. 1213(1 for the pseudocritical coupling (Eq. (HJ 
and the specific heat through Xij i Xee and Xie 
(Eq. Ell as shown in Fig. I2I3I Result for the crit- 
ical exponents are summarized in Table |21 These 
estimates tends to be greater than the expected 
values for 0(4) {yt = 1.34, yt = 2.49), 0(2) 
{yt = 1.49, yh = 2.48) and mean field {yt = 3/2, 
yh = 9/4). Moreover the exponent from p seems 
to indicate a first order phase transition. As a 
consistency check one can fit the peak of Xm to 
read back from data the assumed magnetic ex- 
ponent yh (Fig. El . Also in this case we note a 
deviation from the starting value yh = 2.49. 

We tested the 0(4) scaling behavior and found 
it inconsistent with our data. A first order tran- 
sition cannot be excluded but our initial assump- 
tion for yh flaws the analysis. Moreover the 
staggered fermion action at A^t = 4 is probably 
plagued by large lattice artifacts that suppress 
the critical behavior. 
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Figure 3. External lines show the 1-a confidence 
region for yh 
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